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Introduction: SI Units 

m = meter, s = second, kg = kilogram, N = Newton, J = Joule 

. Distance x [ m ], Time t [ s ] 

. Speed (velocity) v = ^ or ir = ^ 

v [ m per s = m/s = m s _1 ] 

. Acceleration a =% = ^ or a = ir 

a [ m/s 2 = m s~ 2 ] 

. Force F = ma 

m [ kg ], F [ N = kg m s~ 2 ] 

. Work W = J F dx or w = fax 

W [ J = N m ] 

. Enerqy W = AE = r F dx (Kinetic Energy, Potential Energy, Mass Energy) 

E [J] or [eV] [ eV = 1 .602 x 1(D 19 J ] 
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Introduction: Areal density 


Ordinary thickness doesn’t tell you amount of material traversed 

- d = 20 cm Water vs. d = 20 cm Aluminum 

more material traversed in Al, even though distance d same 

Areal density x = dp 

- x [g/cm 2 ] 

- d [cm] 

-p [g/cm 3 ] 

x = 20 g/cm 2 

- Pwater = 1 g/cm 3 p Ai = 2.7 g/cm 3 

- c/water = 20 cm d A i = 7.4 cm 

Same amount of material traversed in 20 cm Water and 7.4 cm Al 

•^ISS ~ 15 g/cirr, X\|ars Atmosph ~ 20 g/cm - , X| -(I'll \inioxph ~ 1000 g C ] 1 1 
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Introduction: Relativity 


E = me 2 

E = T + m 0 c 2 
E = Total energy 

- T = Kinetic energy, often written E !!!!!!!!! 

- m 0 = Rest mass energy 

- m = 7A77 0 , =7 E = 7A77 0 C 2 

-'y=— r J==, /3 = V/C => 7 = -J= 

1 vmuc ) 2 ' sf ^ W - 

Note: T = E - m 0 c 2 / \mv 2 
k, ±mv 2 if v < c 


John Norbury (NASA Langley) 


Space Radiation Transport 


Friday June 19, 2015 


5/72 


Introduction: MeV / nucleon 


Kinetic energy: people write E, but they mean T ! ! ! 

E = T + m 0 c 2 = me 2 = ~/m 0 c 2 

Calculate speed 

1 = ~ = -J— + i 

N /1-(|//C) 2 ' m o c 2 

t [MeV/n] => T = AT 

m Q = Am n , m n = m nacle0B = 938 MeV/c 2 

7 = /no^ + 1 = A^ + J[ = i + 1 

E.g. What is speed of 1000 MeV/n 56 Fe projectile? 
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Cross Section 


o Typical scattering experiment: beam of particles fired into a target 
9 Surrounding detectors detect the particles coming out of reaction 



[Lederman & Schramm, From quarks to the cosmos, Freeman, New York, 1989] 


John Norbury (NASA Langley) 


Space Radiation Transport 


Friday June 19, 2015 


7/72 


Cross Section 


9 Let’s say that some protons are detected 

9 If one just counts the total number N of protons emitted per unit 
time then this is related to the total cross section by 


N = Ca 


9 C = luminosity of beam 

9 a = cross section - related to probability that reaction will occur 
(that a proton will be emitted) 

9 Above formula: more particles will be detected if reaction 
probability is large or if we have very intense beam 
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Cross Section 


o Luminosity typically measured in cm~ 2 s _1 
® E.g. Fermilab proton beam (2000) was about C = 10 32 cmr 2 s _1 
® Cross section is measured in units of area namely cm 2 
® N = Co gives units of counts per second 
® Note that 

barn = 10~ 28 m 2 = 10~ 24 cm 2 

and a picobarn (pb) is therefore (pico = 10~ 12 ) 

pb = 10~ 36 cm 2 

® Fermilab luminosity can be re-written 

£ = 10 32 cm~ 2 s _1 = 10 -4 pb -1 s -1 
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Cross Section 


® Sometimes accelerator administrators use integrated luminosity 
Am, which is simply luminosity multiplied by time 


® E.g. if Fermilab ran continuously over a year, this would 
correspond to an integrated luminosity of 

AnT' 1 = 10 -4 pb -1 s -1 x 365 days « 3000pb _1 

® If the accelerator only delivered say 1500 pb _1 then that would 
mean it was shut down for half the time 
® Integrated luminosity related to integrated number of counts N int 
(total number of counts, not number of counts per unit time) by 

N int = Nt 

where N is number of counts per unit time from before N = Ca 
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Cross Section 


® Ice puck collisions: 

- Kinetic energy conserved if surface frictionless 

® Particle collisions 

- Often kinetic energy not conserved 

- Especially if new particles are created, e.g. p+p -> p+p+7r° 

- Or particles are lost (absorbed), e.g. p + 12 C -a 13 N 

® Total energy E always conserved (just like momentum & charge) 

® Elastic collision = Kinetic energy conserved 

- Does not mean that incident KE remains same 

- E.g. Neutron can undergo elastic collision => impart KE to proton 

® Inelastic collision = Kinetic energy not conserved 
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Cross Section 


9 Total CrOSS section - second definition of the word “total” 

r^total — ^elastic "T ^inelastic "T ^absorption 

9 Scattering cross section is defined as 

^scattering — ^elastic "T ^inelastic 
= ^' ^total — ^scattering "T ^absorption 

9 Reaction cross section is defined as 


^reaction — ^inelastic "T ^absorption 
^ ^total — ^elastic "T ^reaction 

9 Some authors write (including me!) 

^total — ^"elastic "T ^inelastic 
Where <7 reac ti 0 n = ^inelastic 
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Cross Section 

[Lamarsh, Nuclear reactor theory, Addison-Wesley, Reading, Massachusetts, 1966] 

. Elastic scattering =7 Kinetic Energy (KE) conserved 

- Nucleus unchanged in isotopic composition or internal energy 

- KE of neutron and nucleus can change (but total KE conserved) 

- Neutron can transfer KE to nucleus 

- Neutron reappears after interaction (n,n) 

. Inelastic scattering => KE not conserved 

- Nucleus left in excited state 

- Neutron reappears after interaction (n,n') 

Neutron scattered elastically (n,n) or inelastically (n,n') 

- Neutron reappears after interaction 

- Absorption reaction =>■ Neutron disappears 

Examples (n, 7 ), (n,p), (n ,a), fission (by convention) 

Notation: (a,b) means a + target b + anything 
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Cross Section 

[Lamarsh, Nuclear reactor theory, Addison-Wesley, Reading, Massachusetts, 1966] 

- Definitions differ from previous slides, but Lamarsh definitions clearer 

^total = ^elastic(n,n) ^inelastic (n,n') ^( 11 , 7 ) “1“ ^(n,p) “1“ ^(n,2n) ^(n,3n) 

- l _( ^(n,pn) ^(n,a) “1“ ^(n, fission) “ 1 “ ' ' ' 

^total — ^elastic ^nonelastic 


^nonelastic 


^inelastic + ^(n/y) + a (n,p) + a (n,2n) + a (n,3n) 
- l - 0"(n,pn) ^(n,a) ^(n, fission) “1“ ' ' ' 


^absorption 


^(n, 7 ) “1“ ^(n,p) “1“ ^(n,o;) “1“ ^(n, fission) H“ ' ' ' 


^nonelastic — ^inelastic H” ^absorption H - ^(n,2n) ^(n,3n) “1“ ^(n,pn) ' 

Note: C^nonelastic 7 ^ ^inelastic 


John Norbury (NASA Langley) 


Space Radiation Transport 


Friday June 19, 2015 


14/72 


(mb) 


Total Cross Section 


Projectile + Target Detected particles 

Count total number of detected particles 
Total cross section a 

- Easiest to measure, hardest to calculate (from QFT) 

Ne + NaF -> n* + X 

140 
120 
100 
80 
60 
40 
20 
0 

0 0.5 1 1.5 2 2.5 

Tlab (AGeV) 

[Norbury, Nucl. Inst. Meth. B 254, 187, 2007], [Norbury, Nucl. Inst. Meth. B 267, 1209, 2009] 



p + p -> 7T + X 



Plab(GeV) 
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Single Differential Cross Section 


Count particles at different angles (e.g. e = o o ,30°,60°,90°) 

- Angular differential cross section ^ = a(9) (Angular distribution) 

- CT = %d9 = a(9) 69 dn = 27 r sin 9d9 

- Harder to measure 



John Norbury (NASA Langley) 


Space Radiation Transport 


Friday June 19, 2015 


16/72 


Single Differential Cross Section 


Count particles of different energies (e.g. e = 5 Mev, 10 Mev, 20 MeV) 

- Energy differential cross section d f E = a(E ) (Spectral distribution) 
Energy = E = Kinetic Energy!! 

-°=sih% dE =s!r< E )‘ iE 



[Norbury, Physical Review C 37, 407,1988] 
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Double Differential Cross Section 

Count particles of different energies and at different angles 

e.g. ( e = 30° for £ = 5 MeV, 1 0 MeV) and (Q = 60° for E = 5 MeV, 1 0 MeV) 

- Double differential cross section = a(E,9) 

-o - l^dEdfl- j o(E.I))dEd(l 

= -£ = S£m d <> = / 

- Hardest to measure, Easiest to calculate (from QFT) 



[Norbury et al., NASA Technical Paper 2600, 1986] 
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Total Cross Section - Two uses of word “total” 


® Total = Elastic + Inelastic (or Total = Elastic + Inelastic + Absorption) 

Example: Ctotal = ^elastic T <Tinelastic 

Example: = do ^ lasllc + dc7 y ic 

® Total = Non-differential 

«■ = / w df> = / % dE = / & dEde — f a(E,6) dEdO 

Example: <T lnda>llc = / dS 

® cr “total - total” ct t otai = f ^jf- do = f f* + / d(T ™^ 
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Recap: Double Differential Cross Section 


• Often called 3-dimensional (3d) 

• Very precise information 

• Function of both energy & angle 

• Ultimate test of nuclear theory 

• Needed for 3d transport 
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Recap: Single Differential Cross Section 


a(E) = fd9a(E,9) % = / d6 

= tr(E, E') later in BTE 

a(9) = JdEa(E,9) % = f dE 

* Less precise information 

• Function of energy OR angle 


C/ 2 er 

dEde 

d 2 a 

dEde 


John Norbury (NASA Langley) 


Space Radiation Transport 


Friday June 19, 2015 


21/72 


Recap: Total 


Cross Section 


<7 = J dE a(E) = J d9 a (6) = f dEd6 a{E : 6) 

= cr(E) later in BTE 

• Very crude 

• No angle or energy dependence 

• Just total number of counts 
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Cross Section: Exclusive vs. Inclusive 


EXCLUSIVE: A + B^ C+D+E+F 

• Ultimate theoretical test 

• Hard to measure 

INCLUSIVE: A + B^ C + anything 

• Very crude theoretical test 

• Easy to measure 
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Bethe - Bloch equation 


[RMP1260, NCRP 153-62, Jackson628, Alpen369, NASARP1 257-34,40,73, RPP] 
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Boltzmann Transport Equation 

1 -DIMENSIONAL TRANSPORT EQUATION 

A-A s , (e) + 5 ,. (£) 

Solve for 

4 r,{x , E) = differential fluence of isotope / at position x with energy E 
07(E) = cf total = isotopic total cross sect, for isotope / with energy E 

O'total — ^elastic A ^nonelastic — ^elastic A ^inelastic A ^absorption A toother 


/*oo 

Mx, E ) = H j dE ' ^ E - E ') M x ’ E ') 


(jjj(E. E') = inclusive isotopic energy (single) differential cross section 
for producing isotope / with energy E from isotope j with energy E' 

Different materials: a are cross sections for specific target material 

Note: £s,(E)*(x, E) = £[S,(E)0,(x, £)] 
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Boltzmann Transport Equation - Units 


fl -dimensional Transport Equation 


A-A s , (E) + 5i(E) 

M x > e ) = J2 i 

j J 

roo 

1 dE' &ij{E, E') (j)j(x, E') 


® </> ~ [ anything ], x ~ [ g/cm 2 ] => 5-,(£) ~ [ cm 2 /g ] 

« 1 2 g of 12 C contains N A molecules, where N A = 6.02 x 1 0 23 

12 C: number density = number per gram = p [ #/g ] = N A /12 [ g^ 1 ] 
In general, p [ #/g ] = N a /A [ g- 1 ] 


® Thus, 
« S,(E) 


, n = # atoms /cm 3 (number density), p = “normal” mass density [ g/cm 3 ] 
5 ■ = pa ~ [ cm 2 /g ] = [ g _1 cm 2 ] 

= §=> wSm = same units as Note: L = f [MsY], but 


L 

P 


[ MeV^] 
L g/cm z J 


« aij{E,E')= d -^f> ~ [ g' 1 cm 2 E- 1 ] 

® Let differential fluence E) ~ [ # cm -2 E _1 ] 

® E ~ [ anything ], because cancel everywhere, Typically E ~ [ A MeV ] 
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Boltzmann Transport Equation 


3-dimensional Transport Equation 


q v - /e s,(e) + s,(e) 

0,(r,n,E) = £ | ofE'cfn'ffyin.n'.E.E'^/lr.n'.E') 
/ J 


Solve for </>,( r, £2, E ) = differential fluence of isotope / at r moving in 
direction £2 with energy E 

07 (E) = isotopic total absorption cross sect, of isotope / with energy E 

= er before 


<5yy(£2, £2', E. E') = inclusive isotopic double differential cross section for 
producing isotope / moving in direction £2 with energy E from sotope j 
moving in direction £2' with energy E' 

= cr( O, E) before 
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Dose 


® Absorbed dose pcruss] D = ^ 

Units [ Gy = J kg -1 ] 

de = mean energy deposited in mass dm 

9 Dose deposited at position x obtained by ^ over all particles j 


Dose CALCULATED WITH DIFFERENTIAL FLUENCE </> 


9 Sj(E) = stopping power as function of particle kinetic energy E 
9 (j)j(x, E) = differential particle fluence at position x 
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Dose 


® cf)j(x, E ) = differential particle fluence at position x [ cm2 M eV 1 
® <&j = particle fluence (integrated) [ ] (Borak, Nelson lectures) 

good for a beam of fixed Z, fixed E, then LET fixed 


Dose CALCULATED WITH INTEGRAL FLUENCE <t> 

If L has units [ ]: 

L g/cm z J 

D = L<t> 


[— ] A = [—]-[-] 

g/cm 2 cm 2 g kg 


If L has units [ ]: 


D = — 0 

P 


MeV/cm # 
g/cm 3 cm 2 


MeV T 

[ ] ~ [— ] (Borak, Nelson) 

g kg 
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Dose - Unit conversion factor 


Assume L has units [ ^ ] (Borak, Nelson) 

D = -cD 
P 

keV / urn keV cm 

~ — 7 — 3- cm = 

g/ cm j g pm 

10 3 eV 10 _2 m _ 10 eV _ 10 1.602 x 10 19 J 

10 _3 kg 10 _6 m kg kg 

= 1.602 x 10' 9 - 

kg 

= 1.602 x 10~ 9 Gy = 1.602 x 10~ 7 cGy 


Dose conversion factor 

D = k 0 = 1.602 x 10~ 9 Gy with L ~ [ g ] , 0 ~ [cm' 2 ], p ~ [^] 
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Dose Equivalent 


® Stopping power « LET, L = dE/dx 
Units [ keV/^m ] 

9 Dose at a point becomes 


D{x) = '£ dE Lj(E)c/)j(x, E) 


DOSE EQUIVALENT calculated with differential fluence <p 


"(*) = £/ dE Q(Lj) Lj(E) <t>j(x, E) 


[Sv] 


Q(Lj) = quality factor, function of LET, which is function of E 

=> Q(Lj) = Q(Lj(E)) 
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Dose Equivalent 

Previous equation gives dD/dE = Lj(E)(j)j(x, E), so that 

/•OO 

»(*) = £/ dDQ(L , ) 
i Jo 

Defining hcru93, p .5] 

_ dD 
L= ~di 

with units of m kg -1 , gives 

/*oo 

H(x) = J2 dLD L Q(Lj ) 

i Jo 

Note that it is incorrect to write 

roo 

HM/E/, dLD(L) Q(Lj) 
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Effective Dose 


» Different tissues have different radiosensitivity 


Effective Dose 

E = ^w t H t 

[Sv] 

T 
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Summary 


Cross Sections 07 , ay PHYSICS (NUCFRG, QMSFRG) 

Boltzmann TRANSPORT Eqn ^ FLUENCE fa [ #/cm 2 /MeV ] (HZETRN) 

" r\ r\ ~l pOO 

aj - ggS,(e) + cn(E ) J <pj(x. e) = Y, ] e ne ^(e. r> <p,(x, e') 


Dose D(x) 
Dose Equivalent H(x) 

POO 

= ]T / dE Si(E) </>/(x, E) [Gy = J/kg] 

/ - 7o 

POO 

= J2 dE Q{l~i{E)) Li(E) </>;(x, E) [Sv] 

; ^ 0 


Effective Dose F = ' 1/1/7- Hj- [Sv] => Risk 


BIOLOGY Q(Lj),w T 

^ k 5 ^ 
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PlONS 


Leptons spin = 1/2 Quarks spin 

= 1/2 

Flavor 

Mass 

GeV/c 2 

Electric 

charge 

Flavor 

Approx. 

Mass 

GeV/c 2 

Electric 

charge 

K nffio* ( 0 - 0.1 3 ) X1 0-9 

6 electron 0.000511 

0 

-1 

u p 

(\j down 

0.002 

0.005 

2/3 

- 1/3 

y middle 
sjff' neutrino* 

( 0 . 009 - 0 . 13 ) x10“ 9 

0 

CJ charm 

1.3 

2/3 

H muon 

0.106 

-1 

S) strange 

0.1 

- 1/3 

y heaviest 
neutrino* 

( 0 . 04 - 0 . 14 ) x10- 9 

0 

X* top 

173 

2/3 

. X tau 

1.777 

-1 

Jj) bottom 

4.2 

- 1/3 j 


Neutron, Proton = 3 quarks Pion = 2 quarks 

[http://education.web.cern.ch, 201 4] 
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PlONS 


[Slaba, Blattnig, Reddell, Bahadori, Norman, Badavi, Advances Space Research 52, 62, 2013] 



Flowchart of HZE-7 t/EM transport model 


John Norbury (NASA Langley) 


Space Radiation Transport 


Friday June 19, 2015 


36/72 










PlONS MAKE LARGE CONTRIBUTIONS TO DOSE 


[Slaba, Mertens, Blattnig, NASA TP-2013-217983] 

[Norman, Blattnig, De Angelis, Badavi, Norbury: Advances Space Research 50, 146, 2012] 
[Slaba, Blattnig, Reddell, Bahadori, Norman, Badavi: Advances Space Research 52, 62, 2013] 
[Aghara, Blattnig, Norbury, Singleterry: Nucl. Inst. Meth. B 267, 1115, 2009] 



1 Pions + Muons + EM 
Ions + Neutrons 


MDU=Mobile Dosimetry Unit (ISS) 

Important discovery in space radiation 

John Norbury (NASA Langley) Space Radiation Transport Friday June 19, 2015 


37/72 


Minimum in Dose Equivalent vs. Depth 

[Blattnig, Slaba, Bahadori, Norman, Clowdsley, Space Radiation Investigators’ Workshop, Galveston, TX, 2014] 

Previous design paradigm 

Exposures from GCR show little variation past 40 g/cm 2 

- Increased shielding slightly decreases exposures 

- Transport performed with HZETRN (straight ahead transport with no pions) 




BON2010 used to generate Oct 1976 GCR environment. 
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Minimum in Dose Equivalent vs. Depth 



Figure 3. Dose equivalent (left pane) and effective dose (right pane) as a function of shield thickness on the 

Martian surface. 
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Minimum in Dose Equivalent vs. Depth 




Figure 4. Comparison of surface exposure to !4 free space exposure as a function of shield thickness. 
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Minimum in Dose Equivalent vs. Depth 

[Blattnig, Slaba, Bahadori, Norman, Clowdsley, Space Radiation Investigators' Workshop, Galveston, TX, 2014] 

New design paradigm 

If forward/backward (FB) neutron transport and pions turned ON 

- Minimum in dose equivalent vs. depth near 40 g/cm 2 

- Increased shielding increases exposure 

- Material optimization more important than previously thought 




BON2010 used to generate Oct 1976 GCR environment. 
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Minimum in Dose Equivalent vs. Depth 


[Blattnig, Slaba, Bahadori, Norman, Clowdsley, Space Radiation Investigators’ Workshop, Galveston, TX, 2014] 

Dose equiv in Al is result of neutron production & back-scatter 
in shielding behind target point 

Back scattered neutrons produce high LET target fragments 
with short ranges 


Front shield Back shield 



t 


Target point 
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Minimum in Dose Equivalent vs. Depth 


[Blattnig, Slaba, Bahadori, Norman, Clowdsley, Space Radiation Investigators’ Workshop, Galveston, TX, 2014] 

Effective dose - body self shielding moderates neutron field 
& reduces impact of backward component 
FB neutron transport accounts for multiple elastic collisions (moderation) 


Front shield/tissue Back shield/tissue 
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Minimum in Dose Equivalent vs. Depth 


[Blattnig, Slaba, Bahadori, Norman, Clowdsley, Space Radiation Investigators’ Workshop, Galveston, TX, 2014] 

Pion/EM contribution to dose equivalent smaller than contribution from 
neutrons, but still noticeable 

- Neutron contribution mainly from high LET target fragments (avgQ~20) 

Pion/EM avg Q ~ 1 lOQQr c l-ji- uvtttdat.cd , , m,,. 

Solid lines: HZETRN (FB neutrons + tc/EM) 

Dash-dot lines: HZETRN (FB neutrons) 

Dotted lines: HZETRN (Straight ahead) 



20 40 60 80 

Shield thickness (g/cnT) 


100 


John Norbury (NASA Langley) 


Space Radiation Transport 


Friday June 19, 2015 


44/72 


Minimum in Dose Equivalent vs. Depth 

[Blattnig, Slaba, Bahadori, Norman, Clowdsley, Space Radiation Investigators’ Workshop, Galveston, TX, 2014] 

Pion/EM contribution to effective dose more noticeable because 
tissue moderates/attenuates contribution from neutron field 
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Minimum in Dose Equivalent vs. Depth 

[Blattnig, Slaba, Bahadori, Norman, Clowdsley, Space Radiation Investigators’ Workshop, Galveston, TX, 2014] 

Verification 

Preliminary comparison to Geant4 


Front Al BackAI 

shield (x) shield (x) 




BON2010 used to generate solar minimum GCR 
environment (<{> = 475 MV). 
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Minimum in Dose Equivalent vs. Depth 


[Blattnig, Slaba, Bahadori, Norman, Clowdsley, Space Radiation Investigators' Workshop, Galveston, TX, 2014] 

Summary & Conclusions 

• New improvements to transport codes indicate potentially significantly 
different exposure vs. depth behavior in the region of 30-100 g/cm 2 

• Significant work has been done to confirm these findings but more is 
needed including: 

- More extensive benchmarking with Monte Carlo codes 

- Thick target experiments 

• Vehicle and system designs may need to adapt to these new results 
including: 

- Renewed emphasis on materials optimization 

- Changes in configurations to avoid thicknesses where exposures are maximal 
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Internal Environment 


NASA space radiation has focused on heavy ions (Fe) 

- Biology: majority of experiments at NSRL 

- Nuclear Physics: QMSFRG based on Eikonal approx 

- Transport: Id HZETRN based on Straight-Ahead approx 


- Reason 
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Relative Abundance 


Internal Environment 


(FREE SPACE) 



Atomic Number (Z) 


[NCRP Report No. 98] 
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Contribution (% ) 


Internal Environment 


(FREE SPACE) 



GCR environment. Fluence, Dose, Dose equivalent. 
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Internal Environment 


(5g/cm 2 ) 


c 

o 

• 4 — ' 

-Q 


C 

o 


O 



0 5 10 15 20 25 30 


Charge Number 

Relative contribution in fluence, dose and dose equivalent of different elements in the GCR 
spectrum. Calculation is an average over 1 year in solar minimum behind 5 g/cm 2 Al shielding. 

[Durante & Cucinotta, Rev. Mod. Phys. 83, 1245, 2011] 
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Internal Environment 


BUT we don’t have naked astronauts in free space . . . 

... and usually not thinly shielded 
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Internal Environment 


[Walker, Townsend, Norbury, Advances in Space Research 51, 1792, 2013] 



■ Z=0 1Z=1 ■ Z=2 1 3<Z<10 ■ 11<Z<20 ■ 21<Z<28 

Percent contribution to BFO dose equivalent by charge group. 

Light ions & neutrons dominate! 
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Internal Environment 


NASA space radiation has focused on heavy ions (Fe) 

- OK for naked astronaut or thin shielding 

- Biology: more neutron, light ion experiments at NSRL 

- Not saying Fleavy ions not important: need both heavy + light 

Neutrons & Light ions deflected at large angles 

- Nuclear Physics: QMSFRG based on Eikonal approx 

- Transport: Id HZETRN based on Straight-Ahead approx 

- Both not adequate for neutrons, light ions 

Langley: 3d nuclear physics & 3d transport for neutrons, light ions 

Neutrons large Q, light ions highly penetrating 

- Accurate description ->• lower risk 
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GCR Sensitivity Analysis: external environment 


Heavy Ions in external environment 

Produce the light ions and neutrons in the internal environment 
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GCR Sensitivity Analysis 


ACE = Advanced Composition Explorer 
Earth-Sun LI (1.5 million km from Earth) 

- CRIS Cosmic Ray Isotope Spectrometer 

AMS = Alpha Magnetic Spectrometer (on board ISS) 

- Search for dark matter, antimatter etc. 

- Also cosmic rays 

NSRL = NASA Space Radiation Lab 

Brookhaven National Lab on Long Island 

- p - Fe 50 MeV/n - 1 GeV/n 
-uptoAu 165 MeV/n 

FAIR = Facility for Antiproton and Ion Research 

- GSI Gesellschaft fur Schwerionenforschung (Darmstadt) 

- Ne 45 GeV/n, U 35 GeV/n 
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Differential effective dose (mSv - MeV/n‘-ycar) 


GCR Sensitivity Analysis - 20 g/cm 2 


[Slaba & Blattnig, Space Weather, in press, 2014] 




Figure 2. Differential effective dose rate as a function of boundary kinetic energy behind 20 g/cm 2 of aluminum 
exposed to solar minimum conditions described by BON2010 model. Results for specific ions have been scaled to 
improve plot clarity. The location of the peak distribution values along the horizontal axis indicates which boundary 
energies are most important to effective dose behind shielding. 

GCR spectrum 90% effective dose > 500 MeV/n (ACE not helpful) 

Friday June 19, 2015 


John Norbury (NASA Langley) 


Space Radiation Transport 


57/72 


GCR Sensitivity Analysis - 0 g/cm 2 


Energy Range 


Satellite Accelerator 


Ei < 250 MeV/n ACE, AMS 

E 2 250 - 500 MeV/n ACE, AMS 

E 3 500- ^ 500 MeV/n AMS 

E 4 1500 -4000 Me V/n AMS 

E 5 > 4000 MeV/n AMS 


NSRL 

NSRL 

NSRL, FAIR 
FAIR 
FAIR 


[Slaba & Blattnig, Space Weather, in press, 2014] 

Table 1. Relative contribution (xlOO) of GCR boundary energy and charge groups to effective dose with no 
shielding. A value of 0.0 indicates that the relative contribution is less than 0.1%. The BON2010 GCR model was 
used for these results d uring solar minimum conditions. 



a 

E 2 

£3 

®4 

A 

Total 

Z= 1 

2.1 

3.7 

11.6 

10.5 

7.7 

35.6 

Z = 2 

2.1 

1.4 

2.5 

1.6 

0.9 

8.4 

Z= 3-10 

4.1 

5.5 

2.6 

1.0 

0.6 

13.8 

Z= 11-20 

1.4 

6.0 

7.7 

2.7 

1.5 

19.3 

Z = 21-28 

0.4 

3.0 

11.1 

5.5 

3.0 

22.9 

Totals 

10.0 

19.5 

35.6 

21.3 

13.6 

100.0 

e 3 + e 4 

+ e 5 

= 70% 

e 4 + E5 = 

35% 
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GCR Sensitivity Analysis - 20 g/cm 2 


Energy 

Range 

Satellite 

Accelerator 

E, 

< 250 MeV/n 

ACE, AMS 

NSRL 

e 2 

250- 500 MeV/n 

ACE, AMS 

NSRL 

e 3 

500- 1500 MeV/n 

AMS 

NSRL, FAIR 

E, 

1 500 - 4000 Me V /n 

AMS 

FAIR 

e 5 

> 4000 MeV/n 

AMS 

FAIR 


[Slaba & Blattnig, Space Weather, in press, 2014] 

Table 2. Relative contribution (xlOO) of GCR boundary energy and charge groups to effective dose with 20 g/cm 2 
aluminum shielding. A value of 0.0 indicates that the relative contribution is less than 0.1%. The BON2010 GCR 
model was used for the se results during solar minimum conditions. 



E i 


E 3 

a 

£5 

Total 

Z = 1 

1.2 

5.4 

18.2 

18.4 

14.8 

58.1 

Z = 2 

1.2 

2.2 

4.1 

2.9 

1.7 

12.2 

Z= 3-10 

0.0 

3.3 

3.8 

1.3 

0.8 

9.1 

Z = 11-20 

0.0 

0.2 

6.6 

2.0 

1.1 

10.0 

Z= 21-28 

0.0 

0.0 

4.7 

3.8 

2.1 

10.6 

Totals 

2.5 

11.1 

37.4 

28.4 

20.5 

100.0 

E 3 + E 4 + E 5 

= 86 % 

E 4 + E 5 = 

49% 
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GCR Sensitivity Analysis - 40 g/cm 2 


Energy 

Range 

Satellite 

Accelerator 

E, 

< 250 MeV/n 

ACE, AMS 

NSRL 

e 2 

250- 500 MeV/n 

ACE, AMS 

NSRL 

e 3 

500- 1500 MeV/n 

AMS 

NSRL, FAIR 

E, 

1500 -4000 Me V/n 

AMS 

FAIR 

e 5 

> 4000 MeV/n 

AMS 

FAIR 


[Slaba & Blattnig, Space Weather, in press, 2014] 

Table 3. Relative contribution (xlOO) of GCR boundary energy and charge groups to effective dose with 40 g/cm 2 
aluminum shielding. A value of 0.0 indicates that the relative contribution is less than 0.1%. The BON2010 GCR 
model was used for the se results during solar minimum conditions. 



A 

E 2 

£3 

A 

A 

Total 

Z = 1 

0.2 

5.4 

20.9 

23.0 

20.0 

69.6 

Z = 2 

0.1 

2.5 

4.9 

3.9 

2.5 

14.0 

Z= 3-10 

0.0 

0.5 

3.9 

1.4 

0.9 

6.7 

Z= 11-20 

0.0 

0.0 

2.9 

1.4 

0.8 

5.1 

Z = 21-28 

0.0 

0.0 

1.1 

2.3 

1.2 

4.6 

Totals 

0.3 

8.5 

33.8 

32.0 

25.4 

100.0 

E3 + E4 

+ e 5 

= 91% 

E4 + E5 = 

57% 
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Flux (m 2 s sr MeV u ' 1 


GCR Sensitivity Analysis 
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Flux (m 2 s sr MeV u ' 1 


GCR Sensitivity Analysis 
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Flux (m 2 s sr MeV u ' 1 


GCR Sensitivity Analysis 
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Flux (m 2 s sr MeV u ' 1 


GCR Sensitivity Analysis 
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Flux (m 2 s sr MeV u ' 1 


GCR Sensitivity Analysis 
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Flux (m 2 s sr MeV u ' 1 


GCR Sensitivity Analysis 
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Flux (m 2 s sr MeV u ' 1 


GCR Sensitivity Analysis 
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Flux (m 2 s sr MeV u ' 1 


GCR Sensitivity Analysis 
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GCR Sensitivity Analysis - Rough Conclusions 


No shielding (naked astronaut) 

Satellites: ACE helpful, AMS much better 

Accelerators: NSRL very useful, FAIR very useful 

Moderate shielding (20 g/cm 2 ) 

Satellites: ACE inadequate, AMS very useful 

Accelerators: NSRL useful for low Z, FAIR useful for high Z 

Thicker shielding (40 g/cm 2 ) 

Satellites: ACE useless, AMS very useful 

Accelerators: NSRL inadequate, FAIR very useful 
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GCR Sensitivity Analysis - Overstated Conclusions 


NASA Space Radiation Lab (NSRL): 

Wrong Ions (Z)? 

Internal spectrum neutrons & light ions 
Wrong Flux (dose rate)? 

Flormesis, Bystander, Fractionation, Countermeasures 
Wrong Energy (E)? 

- 60% effective dose > 1 .5 GeV/n (4o g /cm 2 ) 


John Norbury (NASA Langley) 


Space Radiation Transport 


Friday June 19, 2015 


70/72 


Summary 


E = 1 GeV/n, really means T = 1 GeV/n 
Cross Section 

■ Total, Differential, Angular, Spectral, Double Differential 
Transport 

Flux, Fluence, Dose, Dose Equivalent 
Pions make significant contributions to dose 
Possible minimum in Dose equivalent vs. Depth curve 
ACE not sufficient - AMS needed 
NSRL not sufficient - FAIR needed 
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THE END 


john.w.norbury@nasa.gov 
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